Corrections to Scaling in the Droplet Picture of Spin Glasses 
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The energy of a droplet of linear extent I in the "droplet theory" of spin glasses goes as I 6 for large 
I. It is argued that this formula needs to be modified by the addition of a scaling correction 
in order to accurately describe droplet energies at the length scales currently probed in numerical 
simulations. With this simple modification all equilibrium numerical data on Ising spin glasses in 
two, three and four dimensions becomes compatible with the droplet model. 

PACS numbers: 75.50.Lk, 02.60.Pn, 75.40Mg, 75.10.Nr 
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The controversy as to the nature of the ordered state 
of spin glasses still remains to be resolved. There are two 
theories: the "droplet theory" pi E| and the replica 
symmetry breaking (RSB) theory of Parisi 0, || . Many 
numerical studies have been done in an attempt to re- 
solve the controversy with results which are often un- 
comfortable for proponents of both theories. The most 
natural summary of the current numerical situation is 
the TNT {trivial, non-trivial) picture of Refs. 0. On 
this picture, which we shall explain in more detail be- 
low, the Parisi spin overlap function P(q, L) for a system 
of linear dimension L seems to behave as expected from 
the RSB theory i.e. in a non-trivial fashion, while the 
behavior of the link overlap seems to be more in accor- 
dance with droplet model ideas as its variance seems to 
be tending to zero as L increases to infinity (which is 
regarded as trivial). In this paper we shall show that a 
simple and natural correction to scaling term added to 
the usual droplet energy expression can explain the origin 
of the TNT picture. 

In the RSB theory there are low-energy excitations in 
which a finite fraction of the spins in the system are re- 
versed but which only cost a finite amount of energy in 
the thermodynamic limit. The surface of these excita- 
tions is space-filling so that the fractal dimension of their 
surface d s is the same as the space dimension d B. In 
the droplet theory, the lowest energy excitations of linear 
spatial extent I typically costs of order I , where 9 is a 
positive exponent, of order 0.20 in three dimensions 
and 0.70 in four dimensions G3| according to studies of 
the energy to create a domain wall right across the sys- 
tem. Thus in the thermodynamic limit the excitations 
which flip a finite fraction of the spins cost an infinite 
amount of energy and also the surface of these excita- 
tions is not space filling, as d s < d. 

The Hamiltonian which is usually studied in numerical 
work is the Ising spin glass model: 



H — ^ JijSiSj 



(1) 



<ho> 



neighbor sites on the lattice. Typically for equilibrium 
studies at low temperatures (done at temperatures T 
much less than the critical temperature T c in order to di- 
minish the complications from critical point effects p| ), 
L < 14 for 3D and L < 5 for 4D. By studying two real 
replicas of the system, with spins Si and 5 f one can derive 
the Parisi spin overlap function P{q, L), where 



q 



N ^ SlS i 



and study the properties of link overlaps 
Ql 



— — SiSjS^Sj, 

<ij> 



(2) 
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where the sites i lie on the sites of a cubic lattice with 
N = L d sites, Si = ±1, the Jy have a Gaussian distribu- 
tion of zero mean and unit variance and couple nearest- 



where Nf, is the number of mearest-neighbor bonds on 
the lattice. On the droplet picture P(q, L) should tend 
as L — * oo to a delta function at q — qEA, and P(q, L) ~ 

T/L e for all q ^ qEA where the expectation is that 9 = 
9. In fact it is found that 9 w and that P(q, L) looks 
very similar to what is predicted from the mean-field RSB 
picture § in that P(q,L) ^ for all q ^ qEA It is this 
result which advocates of the droplet model have found 
hardest to understand. 

Studies of link overlaps are more in accordance with 
droplet model expectations. Their variance is predicted 
on the droplet model to decrease with L as T/L^ 1 with 
IH = 2(d — d s ) + 9 The most recent investigations 

of link overlaps have been performed by adding a bulk 
perturbation —eqi to 7i and studying its effect in the 
ground state (see below). Using this method Palassini 
and Young @ report that in 3D d s = 2.58 ± 0.02 and 
6' = 0.02 ± 0.03 while in 4D d s = 3.77 ± 0.05 and 9' = 
0.03 ± 0.05. The droplet model expectation is again that 
9 =9, which is clearly not satisfied at the L values 
which can be currently studied. On the RSB picture the 
variance of qi should be non-zero for large L, which is not 
in accord with the results of JtJ. (However, Marinari and 
Parisi jl3| did obtain a finite variance of the link overlap 
in the large L limit by extrapolating it in an ad hoc way 
as 1/L). 

Thus if one accepts the TNT version of the numerical 
results but believes in the droplet model, the task is to 
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understand why 9 ^ 9 at least at the small values of 
L which can be currently reached in equilibrium numer- 
ical studies. We will assume that the droplet picture is 
essentially correct and explain why corrections to scal- 
ing result in 9' appearing to be close to zero for small L 
values. 

It is important to realise what the low-energy excita- 
tions look like in spin glasses. Consider a domain wall 
crossing a system M X M X L. The energy cost of this 
domain on the droplet picture will be of order M and the 
domain wall will be fractal with dimension d s and with 
an area of order M ds . In three or more dimensions the 
interface may have holes throught it. Because it is fractal 
the extent of this wandering is of order M. The wander- 
ing of the interface by an amount of order M affects the 
determination of 9, see Ref. Oj, where it was found that 
the best results were obtained when L>Mas then the 
interface is not affected by its interactions with the ends 
of the system. Suppose now we have two domain walls 
across the system. Then if their separation is large com- 
pared to M they will be unaffected by the presence of 
the other. However, if they are closer together then they 
interfere with each other and their overall energy will be 
greater than if the other one were absent. In other words, 
domain walls effectively repel each other. If they have a 
separation of order I, then the repulsive energy between 
them would be expected to vary with their separation as 

a power law, . No investigations of u> seem to exist 
in the literature. 

Consider now a droplet of linear extent /. It too will 
have a fractal surface described by d s . Pictures of large 
droplets have now appeared (^] , and a systematic investi- 
gation of them is in JT^ | . Their fractal nature ensures that 
in three or more dimensions that they have holes through 
them, giving them a sponge-like appearance. However, 
because the surface of the droplet may wander by a dis- 
tance of order I the energy of the droplet E will be modi- 
fied by its wandering and "collisions" with itself to a form 
which we suppose by analogy with the above is 



E = Al e + Bl~ 



(4) 



The term £?Z~" is a scaling correction to the form of the 
droplet energy at large I and A and B are positive con- 
stants. We shall argue that with this scaling correction 
it is possible to understand the TNT description of the 
numerical scene. The exponent u is only well-defined as 
a correction to scaling exponent when I — > oo, but we 
shall assume that Eq. (||) has utility outside this limit. 
We do not know whether uj = u>, but fortunately there 
already exists numerical data on the value of u> in 3D 
from the work of Lamarcq et al. fl6|| . 

These authors first computed the ground state for both 
N = 6 3 and N — 10 3 . They then chose an arbitrary ref- 
erence spin and flipped it along with a cluster containing 
v — 1 other spins connected to it. They next minimized 



the energy of this cluster by exchange Monte Carlo, but 
with the constraint that the reference spin is held fixed 
and the cluster was always connected and of size v. They 
found the largest extension (mean end-to-end distance) I 
of the cluster and found that for v < 33 its energy E(l) 
varied as where the exponent u> was 0.13±0.02. For 
larger v values they found that E(l) was increasing rather 
than decreasing with I but they were unsure whether this 
might not be an artifact of insufficient numbers of Monte 
Carlo steps. Note that according to Eq. (0), when both 
9 and u> are small E(l) will have a shallow minimum. 
The clusters generated by their procedure were not com- 
pact (that is, v ~ I ), but on the droplet picture, which 
is a scaling picture associated with the zero-temperature 
fixed point l], it is only required that large clusters be 
compact. 

Eq. (0) is more illuminating when expressed in vari- 
ables associated with its minimum at I = Iq: 



Eil) = C 



f(i/i y (i /iy 

\ 8 oj 



(5) 



where C is a positive constant. The position of the min- 
imum is essentially unknown but there is a hint from the 
upturn in E(l) seen in Ref. that Iq might lie between 
7 and 10 lattice spacings in 3D. Then because the values 
of both 9 and ui are small there is only a weak dependence 
of E(l) on I in the region accessible to numerical studies 
(say 14 lattice spacings). The minimum is very shallow 
so that the apparent value of 9, which is 9 , will be close 
to zero. For example, the ratio E(14)/E(7) = 1.006 if l 
is 7, whereas in the absence of the correction to scaling 
term this ratio is 1.15. In 4D the values of ui and Iq are 
undetermined at present, but provided Iq happens to lie 
in the region open to numerically studies, then again 9 
would be small. 

While it is inevitable that exponents in 4D will be dif- 
ficult to determine in numerical studies this is not so for 
2D. There seems to be indirect evidence which supports 
an E(l) as in Eq. ([|) also for 2D systems. In two dimen- 
sions 9 is very accurately determined as systems of 480 2 
can be studied JTtJ . The exponent associated with a sin- 
gle domain wall, that is 9, is -0.282(2) How ever, 
studies where, say, 9 is determined from the effects of 
thermally excited droplets (such as Monte Carlo simula- 
tions of the spin-glass susceptibility) yield an apparent 9 
close to —0.47 |Q . This discrepancy has long been a puz- 
zle, and has prompted suggestions that perhaps different 
exponents describe domain wall energies and droplet en- 
ergies [ fl8| ]. However, Eq. ([I|) with u) ~ 0.47 would seem 
to offer another way of resolving the descrepancy. For 
small values of I, droplets would apparently have ener- 
gies decreasing as as their energies are dominated by 
the correction to scaling term, but at large values of I the 
decrease with I will be slower and be as l~^ e \ as expected 
from the conventional droplet approach. 
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This crossover between the large and small I depen- 
dencies can be seen in the work of Middleton jy^. He 
studied link overlaps as in H and |]l3) , by comparing the 
unperturbed Jy ground state, {S^}, with the e-perturbed 
state, where Jy — > Jy — eSfSj/Nj,. The fraction of link 
values SiSj on nearest-neighbor bonds < ij > which are 
changed by the perturbation is (1 — qi). On the droplet 
model the sample average of this, 1 — qi, is predicted 
to decrease as e/Z w . Using the values 8 = —0.28 and 
d s = 1.27, m k 1.18, but with 8 at its small I value of 
—0.47 the expected value of //; « 0.99. Middleton stud- 
ied the effective value of /ii at length scale L using the 
definition 



// eff 



jKl - gi)] 
d\n(L) ' 



(G) 



At L = 16 jtt; ^ was found to be 0.99 and increased to 
the large I expected value 1.18 when L ~ 200. In other 
words, /i; 6 ^^ seems to interpolate beteen the two limits as 
expected on the basis of our correction to scaling. 

Middleton himself attributed this gradual evolution of 
/ij with L to finite size corrections which can arise 
when studying quantities which have contributions from 
droplets of all sizes up to the system size L. The do- 
main wall energy exponent 8 is free of this problem as 
by definition domain walls only exist on the scale L. He 
showed that even if E(l) varies just as l e for all I, then 
the contributions of droplets on all length scales up to 
L means that approaches its asymptotic value only 
as 1/L These scale averaging corrections decrease 

in 2D with the exponent d — hi ~ 0.82 and are therefore 
unlikely to be important when L ~ 100. In three dimen- 
sions d— fii w 1.3. As in 3D only modest values of L can 
be studied the scale averaging corrections may there be 
a significant effect. 

The effect of the correction to scaling in the formula 
for the droplet energy will have relevance beyond the con- 
fines of numerical simulations. The droplet model pro- 
vides scaling laws for time-dependent quantities such as 
dynamical susceptibilities in terms of a length scale L(t) 
which increases in a logarithmic manner with time t due 
to thermal activation of droplets on the length scale L(t). 
In practice the length scale which can be explored in real 
experiments is rather limited (typically more than 10 lat- 
tice spacings but probably less than 200 |2(J). The fact 
that this length scale is not very large means that correc- 
tions to asymptotic formula will be important and that 
will include the correction to the scaling energy proposed 
in this paper. 

To summarize: there exists a simple correction to the 
usual scaling formula for the energy of a droplet which is 



consistent with the TNT summary of the numerical data 
on spin glasses. It would be useful to have direct studies 
of the two correction to scaling exponents exponents u> 
and u> . In 3D and 4D, information on the value of I at the 
minimum of E(l), Iq, is needed. As perhaps Iq might have 
been already seen in the 3D work of Ref . |l6| , a precise 
determination could be possible with a little extra effort. 
In 2D an accurate determination of both exponents might 
be achievable with the techniques used in Refs. jlTj and 

I would like to thank Alan Bray for many discussions. 
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